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Abstract 

We point out that the permanent confinement in a compact 2 + 1-dimensional U (1) 
Abelian Higgs model is destroyed by matter fields in the fundamental representa- 
tion. The deconfinement transition is Kosterlitz-Thouless like. The dual theory is 
shown to describe a three-dimensional gas of point charges with logarithmic interac- 
tions which arises from an anomalous dimension of the gauge field caused by critical 
matter field fluctuations. The theory is equivalent to a sine- Gordon-like theory in 
2-\-l dimensions with an anomalous gradient energy proportional to Ic". The Callan- 
Symanzik equation is used to demonstrate that this theory has a massless and a 
massive phase. The renormalization group equations for the fugacity y{l) and stiff- 
ness parameter K{1) of the theory show that the renormalization of K{1) induces 
an anomalous scaling dimension r]y of y{l). The stiffness parameter of the theory 
has a universal jump at the transition determined by the dimensionality and rjy. As 
a byproduct of our analysis, we relate the critical coupling of the sine-Gordon-like 
theory to an a priori arbitrary constant that enters into the computation of critical 
exponents in the Abelian Higgs model at the charged infrared-stable fixed point of 
the theory, enabling a determination of this parameter. This facilitates the computa- 
tion of the critical exponent u at the charged fixed point in excellent agreement with 
one-loop renormalization group calculations for the three-dimensional Xy-model, 
thus confirming expectations based on duality transformations. 



Key words: Gauge theories; Confinement; Kosterlitz-Thouless transition 



Email addresses: kleinert@physik.fu-berlin.de (Hagen Kleinert), 
nogueira@physik.fu-berlin.de (Flavio S. Nogueira), 
asle.sudbo@phys.ntnu.no (Asle Sudb0). 



Preprint submitted to Elsevier Science 



1 February 2008 



1 Introduction 



Gauge theories ind = 2+1 dimensions are often considered as effective theories 
of strongly correlated systems in two spatial dimensions at zero temperature 
Prominent examples of systems to which such theories are hoped to be 
applicable are the high-Tc cuprates in the underdoped or undoped regime. In 
the undoped regime it is known that spinor QED3 is an effective low energy 
theory for the quantum Heisenberg antiferromagnet (QHA) 0]. It is hoped 
that one effectively can account for doping by coupling the gauge theory to a 
scalar boson representing the holon part (charge part) of composite Hubbard- 
operators describing projected electrons, which however do not satisfy simple 
fermion commutation relations. Similar effective theories have a long history 
as useful toy- models in high-energy phxsics 0; H; 0] , and have recently been 
suggested to describe neural networks 7]. 

Of particular interest in the physics of strongly correlated systems is the com- 
pact version of the 2 + 1-dimensional Abelian Higgs model with matter fields 
in the fundamental representation. This is the model we shall be concerned 
with in this paper and for which we shall find the results summarized in the 
Abstract. 



1.1 Preliminary remarks 

Our starting point is the following Abelian euclidean field theory of a scalar 
matter field coupled to a massless gauge field 

A = \{d, - tA^M' + ml\M' + y 100^, (1) 

where the subscript zero denotes bare quantities. It corresponds to a theory 
with a Maxwell term 

where = 9^v4° — 9i,y4°, in which the gauge coupling cq goes to infinity. 

This limit implies the constraint = 0, where = (pl d is the boson 
current. 

When deriving effective theories for the t — J model we arrive naturally at a 
compact U{1) lattice gauge field 0. For QHA, the gauge symmetry is larger 
and given by the gauge group 5'?7(2) j^. However, in this reduced U{1) 

formulation is also possible llf. Since this U{1) is a subgroup of SU{2), which 
is a compact group, the f/(l) gauge theory of QHA is necessarily a compact 
Abelian gauge theory. 
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It is well known that a compact U{1) theory of the pure Maxwell type in 
three dimensions confines electric charges permanently jsj. In the literature 
[9] it is also argued that this permanent confinement should be present if an 
additional fermionic field ip coupled to the gauge field by a Lagrangian 

N 

Cf = Y.M^,-^Al)i^.. (3) 

i=l 

This means that the particles represented by the fields ip and 0o never have an 
independent dynamics. In the context of many-body theory, the Dirac fermion 
ip could represent a spinon, while represents a holon. If electric test charges 
were permanently confined in the model, then the spinon and the holon would 
only appear as composite particles. In this case it would be impossible to 
fractionalize the electron, i.e. spin and charge would always remain attached 
to each other. Spin-charge separation is known to occur in 1 -|- 1 dimensions 
[lol |. There fermions can be transmuted into bosons via the so called Jordan- 
Wigner transformation. In 2 + 1 dimensions the situation is less clear, but 
for matter fields in the fundamental representation there is one circumstance 
where spin-charge separation is known rigorously to occur, namely the chiral 



spin liquid state However, the statistics of particles can be changed as 
in 1 + 1 dimensions. In the chiral spin liquid, spinous have anyonic statistics 
described by a Chern-Simons term in the effective gauge theory, which 
reflects the breaking of parity and time reversal symmetry. 

The lack of consensus about spin-charge separation in 2 + 1-dimensional com- 
pact f/(l) matter-coupled gauge theories with matter fields in the fundamental 
representation initiated investigations of other gauge theories for strongly cor- 
related electron systems. One of the most promising candidates is a Z2 gauge 
field coupled to matter fields [3| • Similar ideas leading to electron fractional- 
ization had earlier been presented in the condensed matter literature [l4: isl . 



In 2 + 1 dimensions the Z2 theory has a deconfinement transition j5[. Thus, 
Z2 gauge theories are potentially good candidates for describing spin-charge 
separation without breaking parity and time reversal symmetries. 

The confinement properties of U{1) gauge theories for the cuprates and the 
relation to spin-c harg e separation were recently discussed from various points 
of view [la 12 ; 13 3 • Nayak Q states that in gauge theories of the t — J 



model fermions and bosons interact at infinite (bare) gauge coupling and, for 
this reason, it is necessarily a theory with permanent confinement of slave 
particles. In contrast, Ichinose and Matsui |l^ have argued that the coupling 
to matter fields strongly infiuences the phase structure of the system. In Ref. 
fig^ . it is correctly pointed out that if spin-charge separation occurs, it is not 
necessaryly tied to the notion of confinement- deconfinement of slave particles. 
The picture proposed in Ref. in 2 -|- 1 dimensions is reminiscent of 1 + 1 
dimensions where spinous and holons are solitons and cannot be identified 
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with the slave particles, which are not part of the spectrum • Nagaosa and 
Lee [l7| discuss a compact U{1) gauge theory coupled to bosonic matter field 
in the fundamental representation. They conclude that in = 2 + 1 this theory 
permanently confines electric charges, in contrast to the analysis by Einhorn 
and Savit on the same model 



In a recent letter j20|, we have studied the confining properties of the La- 
grangian (1), as well as the case of a fermionic field ip coupled to a gauge 
field, but with an added Maxwell term. The Lagrangian (1) with a Maxwell 
term corresponds essentially to the model considered by Nagaosa and Lee 
[l7| . though these authors have considered a frozen-amplitude version of the 
model. In Ref. jl^, it was emphasized that an anomalous scaling dimension of 
the gauge field, arising from matter-field fiuctuations, changes the interaction 
between monopoles from 1/r to Inr in three dimensions. It was then argued 
that a monopole-antimonopole unbinding transition similar to the Kosterlitz- 
Thouless (KT) transition takes place, but now in three dimensions. From this, 
we concluded that test charges undergo a deconfinement transition. 

It must be pointed out that the authors of Refs. ^ were looking for a 
transition similar to those encountered in d = 3 + 1, namely ordinary first- or 
second-order phase transitions In Ref. ^^I' ^ duality transformation was 
performed showing that the disorder parameter {(pv) is always different from 
zero, implying that (0) is always zero. This result is essentially correct and is 
perfectly consistent with the scenario in Ref. and explained further in the 
present paper. 



A main result in our letter |20|] is that there exists a non-trivial infrared stable 
fixed point in the theory ind = 2+1 which drives the deconfinement transition. 
There the anomalous dimension of the gauge field is given by r^^ = 1 in 
(i = 2 + 1 j^; 2^1. This result is exact as a consequence of gauge invariance. It 



implies that the non-trivial infrared fixed point arises at an infinite bare gauge 
coupling. To see this, consider the boson-fermion Lagrangian C = Cf + Ci, + 
Cm- Due to gauge invariance, the gauge coupling renormahzes to = Zac^, 
where Za is the wave function renormalization constant of the gauge field. The 
renormalization group (RG) (3 function for the renormalized dimensionless 
gauge coupling a = e^/ fi has the following exact form in 2 + 1 dimensions 



Pa{a, g) = = [7A(a, g) - (4) 

where g is the renormalized dimensionless |0|^ coupling and 7^ = fid In Z^/ dfi. 
Let us assume that there exist non-trivial infrared stable fixed points a* and 
g^,, where the /3 functions /3a and /3g vanish. We have explained in Ref. jl^] why 
such fixed points must exist. (For similar arguments, see Ref. jHI). Moreover, 
large-scale Monte-Carlo simulations have demonstrated explicitly the exis- 
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tence of such a non-trivial fixed point 2^ 2^1 (see also Ref. Esl 'l . Its existence 
has long been assured theoretically by duahty arguments |26l : l27i | (see also 
Section II-B). We shall not repeat the arguments and details here. Instead, we 
focus on the physical consequences of the non-trivial fixed point. 



We would like to stress an important point, pertinent to d = 2 -|- 1 dimensions, 
and quite different from the situation for = 3 + 1. As a — > a^,, the bare 
coupling Cq must tend to infinity. By definition, the above (3 function is given 
at fixed A, ao, and go- Here, A is the ultraviolet cutoff while ao = Cg/A and 
go = uq/A are the dimensionless bare couplings. The fixed point is reached 
for /i — i> 0. Alternatively, the fixed point is reached for A — oo if /i is held 
fixed. However, since ao is fixed it follows that Cg — > oo as A — ^ oo. Thus 



m 



d = 2 + 1, the fixed point theory is at infinite bare gauge coupling. One might 
object that this infinite gauge coupling cannot be relevant for the cuprates 
which have an infinite value of Cq at all scales, not only in the scale invariant 
regime. This is true, but irrelevant as far as the deconfinement transition is 
concerned, which is determined by the non-trivial fixed point structure. The 
situation is analogous to the 0{N) non-linear a model as opposed to the 
0{N) 0^ model. These models are quite different, but agree with each other 
at the critical point jl^; 2^|, thus belonging to the same universality class. In 



our case, the model with the Maxwell term at the fixed point has the same 
correlation functions as the model without it also at the fixed point. 



To summarize the discussion in the above paragraph, the non-compact action 
with no Maxwell term has the same critical behavior as the compact one at 
the critical point corresponding to a non-trivial fixed point, characterized by an 
infinite bare coupling. Had we started from an infinitely weak bare coupling, 
the only fixed point we would have any hope of reaching for d = 2 + 1 would 
be the Gaussian fixed point. 



In Ref. j20| we have pointed out that chiral symmetry breaking can destroy 
the deconfinement in the fermionic case. We want to point out that for the 
combined boson- fermion model, C = Cf + Ci, + Cm, chiral symmetry breaking 
does not spoil the deconfinement transition. Chiral symmetry breaking occurs 
at a lower value of number of fermion flavours Nf, when also bosons are 
present. Kim and Lee js^ claimed that the critical value of Nf is decreased 
by a factor two. Since we have typically ~ 3 and the physical number 
of fermion components in the cuprates is Nf = 2, Kim and Lee argued that 
spin-charge separation would occur at finite doping [j^ . 
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1.2 Anomalous scaling and the potential between test charges 



The high-energy physics hterature is usually concerned with d = A and use 
low-dimensions only in toy models. In condensed matter physics, however, 
2 + 1-dimensional gauge theories are supposed to describe real physical phe- 
nomena such as the anomalous properties of high-Tc superconductors •31'], or 
the physics of QHA ^ |32|]. For d G (2,4] the gauge coupling /5-function may 
be written as 



Pa{oi,g) = ['jA{(y,g) + d- 4:]a. (5) 

Non-trivial fixed points induce an anomalous scaling behavior in the gauge 
field propagator. In the Landau gauge we have that 



D,Ap) = D{p) - , (6) 



with the large distance behavior given by 



The anomalous scaling dimension is given exactly by |21l : |2 

'nA = 7A{a*,g*) = 4:- d. (8) 
Due to the above result, the propagator (7) in configuration space becomes 

Di^) - - (9) 

for all d G (2,4]. The potential between effective electric charges q{R), sepa- 
rated by a large distance R in {d — l)-dimensional space is given by 



VW ~ (10) 



where 



q {R) ~ ~ , (11) 

Va d-A 

and where A is a short distance cutoff. The anomalous scaling in Eq. (11) is 

a consequence of the coupling to matter fields. Due to it, the potential V{R) 

behaves effectively like 1/i? for ci = 3. For ci = 4, it goes like ln(Ai?)/i?, while 
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for (i = 2, it has a confining behavior proportional to R. The regime governed 
by the Gaussian fixed point has q'^{R) = Qo = const, and corresponds to the so 
called Coulomb phase. In this phase, the four- dimensional theory has V{R) = 
Qq/R, whereas V{R) = gglni? for d = 3. We see that the non-trivial infrared 
behavior induces an effective electric potential between test charges similar to 
that which characterizes the Coulomb phase in d = 4. If we extrapolate to 
d = 2, we obtain V{R) = q^R. Note that in = 2, we obtain a confining 
potential irrespective of whether anomalous scaling is taken into account or 
not. 



In compact abelian gauge theories a confined phase is realized by the formation 
of electric flux tubes connecting electric charges. These flux tubes are the dual 



analogs of the magnetic flux tubes connecting magnetic monopoles |33|; |34 



There is a Dirac relation between the effective electric and magnetic charges 



q{R)qm{R) ~ 1- (12) 
Let us consider now the potential between the magnetic charges 



y-^w-'W^-^i^j^- (13) 

From Eq. (11) we see that for c? = 4 the magnetic potential behaves like 
l/[-Rln(Ai?)]. However, for d = 3 we have 



VUR) ~ ^, (14) 

which is self-dual with respect to the potential between electric test charges. 

The Higgs phase for the electric charges corresponds to V{R) ~ const because 
of the gauge field mass gap. The Higgs phase for magnetic test charges, on the 
other hand, is given by V{R) ~ R. In the electric- magnetic duality picture 
jisi li^ l this Higgs phase for magnetic charges is exchanged by the confined 
phase for electric charges. This scenario should be valid for matter fields in 
the adjoint representation. In the absence of matter fields, a compact 2 + 1 
dimensional gauge theory is definitely confined permanently j^. The above 
result shows that if matter fields are present, a deconfined phase is also pos- 
sible. However, if the matter fields are in the fundamental representation, the 
situations controversial EE ; i 0; H 0. Our recent results in Ref. 
seem to be confirmed by the Monte Carlo work in Ref. 0]. The main purpose 
of this paper is to give more details on the scenario proposed in Ref. j20| and 
to describe a theory for a deconfinement transition in Abelian gauge theories 
coupled to matter fields in the fundamental representation. 
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1.3 Outline of the paper 



In Section II, we consider the lattice duality transformations to the 2 + 1- 
dimensional Abelian Higgs lattice (AHL) model, first the non-compact case 
and later the compact case. We then discuss the possible ordinary first- or 
second-order phase transitions these models can have, with matter fields in 
the fundamental representation for the compact case. 

In Section III, we construct the continuum effective Lagrangian and its dual 
counterpart for the compact 2 + 1-dimensional AHL model when matter-fields 
have been integrated out. Because these are central results of the paper, it 
behooves us to announce them here. 



The dual field theory is given by Eq. (51). It represents a description of a 
three dimensional gas of point charges interacting with a logarithmic pair- 
potential, given by Eq. (49). We emphasize that the 3(i In-plasma action of 
Eq. (49) emerges from an underlying matter-coupled gauge theory, Eq. (38), 
by integrating out the fluctuating matter fields and considering the influence 
of critical matter fluctuations on the gauge-field propagator. The result of this 
procedure is the effective theory Eq. (46). Such matter-field fluctuations endow 
the gauge-field propagator with an anomalous scaling dimension 77^ = 4 — d 
(2I : 2^ 1 which in three dimensions alters the interaction between the monopole 
configurations of the gauge-field from a Coulomb-interaction 1/R to a Ini? 
interaction. 



Recall that in contrast to this, in the classic treatment by Polyakov 8] of 
compact three-dimensional QED with no matter fields, the standard three- 
dimensional sine-Gordon field theory with a quadratic gradient term, describ- 
ing the three dimensional Coulomb gas, is obtained. This action is given by, 
in the notation of Eq. (51) 

Ssg = yJ d^^[ipi-d^)ip - 2zocos<^]. (15) 

Polyakov has demonstrated that Eq. (15) has no phase transition, i. e. 
it is always massive. Our Eq. (51) differs drastically from Eq. (15), due the 
presence of an anomalous gradient term. 

In the first part of section IV, we show using the Callan-Symanzik equations, 
that the effective dual Lagrangian Eq. (51) has a massless and a massive phase 
separated at a critical coupling t^- Hence a phase transition must exist. This 
does not by itself suffice to show precisely what sort of phase transition the 
system undergoes, nor does it allow us to construct the correct flow diagram 
of the coupling constants of the problem. It does, however, suffice to show that 
two different phases exist. Since the propagator of the problem is logarithmic 
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in (i = 2 + 1, a Hohenberg-Mermin- Wagner theorem j36| holds. Under such 
circumstances, it is very natural to conjecture that any phase transition in 
the system, if it exists, must be of a topological character. In the second part 
of Section IV, we construct the renormalization group flow equations for the 
problem and show that the phase transition is of a KT-like type. 

In Section V, we consider the connection between the renormalization group 
functions obtained directly from the Abelian Higgs model, and the KT phase 
transition we find in Section IV. The main point here is that we can use the 
value of the critical coupling of the dual effective Lagrangian for the topological 
defects of the gauge field to fix an a priori arbitrary constant which enters into 
evaluating critical exponents for the non-compact Abelian Higgs model. 

In Section VI, we conclude with a summary and outlook. Appendix A dis- 
cusses another type of sine-Gordon theory also exhibiting a KT-like transition 
in three dimensions. In Appendix B, we derive the fiow equations for the stiff- 
ness parameter and the fugacity of the system defined by Eq. (49), and of 
which Eq. (51) is a field theory formulation. In Appendix C, we compute the 
screened effective potential between charges in the insulating phase of the 3d 
In-plasma. In Appendix D, for completeness, we derive the exact equation of 
state for a rf-dimensional In-plasma with no short- distance cutoff and relate 
the singularities in this plasma to the Callan-Symanzik approach of Section 
IV. In Appendix E, we consider, also for completeness, the duality transfor- 
mation of the AHL model with a Chern-Simons term added. This case is of 
interest in the fractional quantum Hall effect ^] and chiral spin liquids 



2 Duality in the abelian Higgs lattice model 



In this section we review the duality approach to the AHL model. Although 
this is a well studied topic ii9; 4,; 2&:-, 2^ , it is worth reviewing it here in 
order to emphazise the differences and similarities between the non-compact 
and compact cases. In particular, we shall discuss the extent to which these 
cases exhibit ordinary first- or second-order phase transitions. The interesting 
case including a Chern-Simons term will be discussed in Appendix D. 

The essential point is that starting from a non-compact or compact AHL 
model, the dual action has the general form 



5'duai = 2 51 hif,M^^{ri - Yj)hi^ - i2'K ^ 1^ • hj, (16) 
where /ij^ G (— oo, oo) and Ij are integer dual link variables. In the non-compact 
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case Ij satisfy the constraint 

V • 1, = 0, 

whereas in the compact case, the right-hand side is nonzero 



(17) 



(18) 



due to monopole charges Qi G Z. The symbol V denotes the gradient vector 
ona simple cubic lattice of unit spacing with components V^/j = — fi. 

2.1 The non-compact case and the "inverted" XY transition 

In the non-compact case, the partition function of the AHL model is given by 



J, /I 



exp(-5). 



where the action S is given by the Villain approximation 



(19) 



Using the identity 



E 



t/2)m^+ia;m 




E 



,(-l/2t)(a;-27rn)2 



following directly from Poisson's formula 

E F{n)^ E / dxF{ 

n=—oo m=—oo 

we obtain 



x)e 



2-nimx 



(20) 



(21) 



(22) 



WdA 



E <^v-m,,o exp \ 

{mi} I i 



+ iKi ■ nij - —iy X Ai) 

^6 



}• 



(23) 



The Kronecker delta in Eq. (23) is generated by the integrations. Now we 
should integrate out the gauge field Aj. The easiest way of performing this 
integration is by the introduction of an auxiliary field hj such that the partition 
function can be rewritten as 
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/oo POQ /"OO 
-GO J —OO J — OO 



Yl dAi^dhi^dk 



y: s{v ■ b,) exp y: 

{Mi} I i 



+ lAi • (bi - V X hi) - —hi + 27riMi • b^ 



(24) 



where a summation by parts has been done to replace hj ■ (V x Aj) by Aj • 
(V X hi), and we have used the Poisson formula (22) to replace the integer 
variables nij by continuum variables bj, at the cost of an additional sum over 
integer variables Mj. We may now integrate out A^ to obtain a delta function 
5{hi — V X hj), after which also bj can be integrated out bj, yielding 



{Mi} 



1 



^{V X hi? 



+ — - 27riMi • (V x h. 



(25) 



Summing the last term in the exponent by parts and going over to integer 
variables Ij = V x Mj, we obtain 



00 • 



n d^ii^ 



(^v i^o exp <^ - ^ 



+ —hi - 2m\i ■ hi 



(26) 



Note that the Kronecker delta constraint above is a direct consequence of our 
change to integer- valued variables. If hj is integrated out we obtain 



Z ^ ZqY <^v-ii,o exp 



-2^(5Yk^.D{vi-v^)l^^ 



(27) 



where the Green function G has the large-distance behavior 

g— \/^e|ri-rj| 



47r|rj - r^'l 



(28) 



The factor Zq in Eq. (27) corresponds to the partition function of a free massive 
gauge boson theory. 

Eq. (27) is the dual representation of the partition function for the non- 
compact AHL model. Due to the constraint V • Ij = 0, the integer links Ij 
form closed loops. 
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By taking the limit e — in Eq. (23), we obtain 



mi , (29) 



which is the loop gas representation of the XY model. If, on the other hand, 
we take the limit (3 —>■ oo in E g. (2 7), we obtain the loop gas representation 
of the "frozen superconductor" 



Z\p=oo = ^^-Ufi exp ( J2 ' (30) 
{U} \ e . / 

which has precisely the same form as in Eq. (29). Therefore, the XY model 
is equivalent to the frozen superconductor, provided the Dirac like relation 
= Air'^f] holds. Eq. (27) is a reformulation of Eq. (19) in terms of the 
topological defects of the model, which are identified as integer-valued vortex 
strings forming closed loops. 

If we consider the phase diagram in the — T-plane (with T = we can 

use Eqs. (29) and (30) to estabhsh the critical points on the axes and T, 
corresponding to T — >■ and — > limits, respectively. From Eq. (29) we see 
that when — > we have a XY critical point on the T-axis. Eq. (30) has 
exactly the same form as Eq. (29), but corresponds to the T — limit. The 
critical point in this limit is therefore = Att'^/Tc, with Tc being the critical 
temperature of the XY transition as described by the Villain approximation. 
This is the so called "inverted" XY transition (IXY) |26| . From the existence 
of these two critical points we can establish a phase diagram where there 
is a critical line connecting them 2^. The ordered superconducting phase 



corresponds to the region < e < ef.. 



2.2 The compact case and the absence of an ordinary phase transition 



In the compact AHL model the gauge field G [— 7r,7r]. The corresponding 
Villain action is now given by 



^ = ^ T.i'^A - - 2^™,^)' + ^ El^M-A V.Aa - 27riV,^)2, (31) 
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and in the partition function we should sum over both integers rii/j, and Ni^. 
Using the identity (21) we obtain 



EE 

{n,} {mi} 



A- 

-j-j- Ll^j^ 



n 



27r 



where 



(32) 



1 

— + im ■ iyOi - Ai) + —Wil + imi ■ (V x A^) 



i 

Now we integrate out and 9i to obtain 



(33) 



E ^V-rii.O^V 
{ni},{mi} 



xmj,n, 



exp 



-E 



2/9 



n, + — m. 



= E exp 

{mi} 

/oo 



exp 



E 



^(Vxh,)' + ^-\il-2m\i-h. 



where from the second to the third hne we used the Poisson formula. Note 
the difference between Eq. (34) and its non-compact counterpart Eq. (26). In 
the latter there is a Kronecker delta constraint V ■ Ij = while in the former 
there is no such a constraint. As we shall see, this difference has important 
consequences. We proceed by integrating out /ij^, thus obtaining the partition 
function 



Zq ^ exp 
{li} 



-2t 



JU 



(35) 



where 



D^,^{ri - rj) = (^S^,^ - "^^j ^i^i - rj), (36) 

{-V^ + Pe^)D{r,-rj)^5ij. (37) 

Due to the constraint V-lj = 0, the term containing V^V;, in Eq. (36) does not 
contribute in the non-compact case, and Eq. (27) results. In the compact case, 
on the other hand, V ■ Ij is completely unconstrained and can take any integer 
value. Thus, in order to bring out the differences and similarities between Eqs. 
(35) and (27), and also to identify the character of the topological defects of 
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Eq. (31) appearing in Eq. (35), we can introduce an auxiliary integer- valued 
scalar field Qi such that V ■ Ij = and rewrite the partition function (35) as 



^0 EE ^vi^Q, exp 

{h}{Q^} 



-27r2/5 ^D(r,-r,) {l^^l^^ + ^Q^Q^ 



. (3J 



Whereas the non-compact theory has only closed vortex lines as topological de- 
fects, the compact case contains also open lines with integer-valued monopoles 
of charge Qi at the ends. 

In the limit /3 — 0, Eq. (38), the vortex loops are frozen out and (38) is the 
dual representation of three-dimensional lattice compact QED Q describing 
a Coulomb gas of monopoles in three dimensions. This is equivalent to a 
sine-Gordon model which is always massive in three dimensions, and leads to 
the well known result that compact QED in three dimensions has permanent 
confinement of electric charges, since the monopole gas will always be in the 
plasma phase. As shown by Polyakov we obtain as a consequence that the 
Wilson loop satisfies the area law. 

As in the non-compact case, the limit ^ corresponds to the Villain form 
of the XY model. Thus, if we consider again a phase diagram in the (e^,T)- 
plane we have that a critical point at exists on the T-axis. However, as we 
shall now show, there is no IXY transition in the compact case. To see this, 
let us take the "frozen" limit /3 — > oo in Eq. (38). The result is 



^ = ^0 E E ^v.i.Q. exp ( E in • (39) 

The sum over Qi is trivially done, Y,{Qi} '^v i^.g^ = 1 after which there is no 
constraint. We are left with a trivial sum over Ij giving Jacobi '(9-f unctions 
■(?3(0,e~^'^ ). Since this function is analytic, there is no phase transition on 
the e^-axis, in contrast to the non-compact case. Thus, at first sight it seems 
that there is no phase transition in the compact AHL model with matter 
fields in the fundamental representation, except for the XF-transition on the 
T-axis. That is, there appears to be no ordinary second- or first-order phase 
transition in the interior of the phase diagram of this model. However, in the 
next Sections we shall derive an effective Lagrangian for the compact Abelian 
Higgs model in 2 + 1 dimensions, which will be shown to nevertheless exhibit 
a topological phase transition of the KT type. 
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3 Effective Lagrangian 



This section is one of tlie central parts of the paper in which we shall derive 
an effective field theory for the compact Abelian Higgs model in = 2 + 1 
dimensions. More precisely, we derive a continuum action, Eq. (51) below, 
for the dual model of the system, obtained after matter fields have been in- 
tegrated out leaving an effective theory for the monopoles of the problem. 
It will turn out that the effective dual lagrangian for the 2 + 1-dimensional 
compact Abelian Higgs model, is described by a theory which has many simi- 
larities to the sine-Gordon theory of Polyakov's pure compact electrodynamics 
in (i = 2 -|- 1 . The crucial difference lies in the fact that the gradient term 
in the dual theory receives an anomalous dimension after the matter-fields 
have been integrated out. It is the presence of this anomalous gradient term 
induced by matter-field fluctuations which eventually will lead to the possi- 
bility of a deconfinement transition in (i = 2 -|- 1, in contrast to the classical 
Polyakov-result of permanent confinement pertaining to the pure gauge theory 



3.1 Three-dimensional compact QED 



Let us consider the Euclidean Maxwell action in three dimensions: 



S = l d'x—Flu. (40) 
where F^^, = d^a^ — d^a^. In order to account for monopoles, we have to 
subtract from F^j, the gauge field of monopoles [s^ 

F*5(x) = 27re^,A(x;L), (41) 

where Sx{x; L) is a delta function on lines L. The dual field strength of = 
efj,uxF^/2 has divergences at the end points of the lines L, say 

d^F^' = 2rm{x) = 2nJ2Q, 6%x - x,), (42) 

i 

where Qi may be arbitrary are integers counting the number of lines ending at 
Xi. The shape of the lines is physically irrelevant. They are the Dirac strings of 
the monopoles at Xj. Under shape deformations, F^l undergoes the monopole 
gauge transformations F^^ —>■ F^l + d^A.^^ — d^A.^ which leave F^^ invariant. 

An ordinary gauge transformation can be used to bring F^^{x) to the form 
F^t = -27re,,xdx J d'y ^^,^_ M y), (43) 
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whose dual field strength is 



F^^^ = -2nd, I d^y^^^^^niy). (AA) 
By substituting F^,^ by F^i, — F^^ in Eq. (40), we obtain the action 

^ = I ^'^i^^A + ^ /^'^/ M3^) 4^|J_^| ^(l/)- (45) 

The action (45) corresponds to the continuum counterpart of the /3 — > limit 
of the lattice action in Eq. (38) describing a Coulomb gas of monopoles. This 
is known to be equivalent to a sine-Gordon action as the one in Eq. (15). In 
three dimensions this theory is always massive and it was shown by Polyakov 
that this implies an area law for the Wilson loop. Thus, electric test charges 
in three-dimensional compact QED are permanently confined. 



3.2 Anomalous three-dimensional compact QED 



When bosonic matter fields are present, the topological defects of the theory 
are vortex loops and vortex lines having monopoles with opposite charges at 
the ends. The vortex lines connecting the monopoles have a line tension a 
which vanishes as the scalar bosons become massless. Thus, when the vortex 
lines become tensionless, we are left with a gas of monopoles. However, the 
anomalous scaling of the gauge field due to matter fields alters the interaction 
between pair of monopoles with respect to the ordinary Coulomb interaction 
case. This will lead us to the anomalous Coulomb gas to be described below. 

From the exact behavior of the critical gauge field propagator we have dis- 
cussed in Section IB, we can write an effective quadratic non-/oca/ Lagrangian 
for the gauge field: 



K ~ 1 

= Y^^'lZ^yil/^^^'' (^^) 

where the constant K = K{a^,,g^) and in the second line of Eq. (46) we have 
rewritten Ca in terms of the dual field strength. Specializing to three dimen- 
sions, we have F, = e,^\F^x/2 and 77^ = 1- After introducing an auxiliary 
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vector field b^, we obtain the equivalent Lagrangian: 



^'a = ^^mV^ K + ih,F,. (47) 

In order to take into account the monopoles, we use the expression for as 
given in Eq. (44). By introducing a new field through = d^ip and using 
integration by parts, we obtain: 



Ca = ^id,v)V^id,^) + ^27rn{xMx). (48) 
Integrating out and using Eq. (42), we obtain the monopole action: 

S^on = 2n^KY,QiQjGix, - x,), (49) 
hi 

where 



Thus, instead of having a standard three-dimensional Coulomb gas with inter- 
action potentials l/\xi — Xj\, we have a three-dimensional gas of point particles 
of charge Qi = ±|(5| (with overall charge- neutrality, see Section IVA) with log- 
arithmic interactions, much akin to the situation one has in two dimensions. 
We emphasize, once more, that this is a result of integrating out matter- 
field fluctuations and considering the effect of critical such fluctuations on the 
gauge-field propagator, which is seen to acquire an anomalous scaling dimen- 
sion from these fluctuations, cf. Eq. (46). It therefore seems plausible, at the 
very least, that one should consider the possibility of having a KT-transition 
of unbinding of monopole-antimonopole pairs, but now in three dimensions. If 
this turns out to be the case, then the confinement-deconfinement transition 
in the 2 + 1-dimensional compact Abelian Higgs model with matter fields in 
the fundamental representation, would be of a to'^ological nature with no local 
order parameter, consistent with previous work 



logic 



We are now ready to state one of the main results of this paper. The system 
defined by Eq. (46) and (49) can be brought into the form of a sine-Gordon 
theory, as in the two-dimensional case, but now with an anomalous propagator, 
whose action is 



^ASG = yJ (fx[ip{-d^f'^ip - 2z^COSip], (51) 

where t = Att'^K and Zq = Att'^K(, with ( being the fugacity of the Coulomb 
gas of monopoles. In Eq. (51), S'asg refers to the action of what we name the 
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anomalous sine-Gordon (ASG) theory, since the cubic power of the propagator 
arises from the anomalous scaling dimension of the gauge field. The manner 
in which the coupling constant t enters in Eq. (51) shows that it regulates the 
stiffness of the phase field ip. Since t (x K, we shall in following sections refer 
to K stiffness parameter. 



Renormalization Group analysis of the anomalous sine-Gordon 
model 



This section is another central part of the paper. Here, we shall consider an 
exact scaling argument applied to Eq. (51). The scaling argument will suf- 
fice to demonstrate that this model has a phase transition. We emphasize 
this as an important point, since recent numerical studies have provided 
strong support for the picture proposed in Ref. jl^ that matter-field coupled 
to compact U{1) gauge fields in c? = 2 -(- 1 lead to a recombination of mag- 
netic monopoles into dipoles. For the dual electric charges, this leads to a 
destruction of permanent confinement, and in Ref. 2^1 it was argued that this 
happened, not through any ordinary first or second order phase transition, 
but rather through a KT-like transition. The authors of Ref. were look- 
ing for more conventional phase transitions, concluding that none were found, 
consistent with the results of Ref. [20]. Having established the existence of a 
phase transition, we then go on to argue that it indeed is of a KT-like type. 
The details are as follows. 



4-1 Callan-Symanzik renormalization group analysis 



Let us consider the renormalization of the anomalous sine-Gordon action de- 
fined by Eq. (51). The infrared divergence is easily studied by considering the 
cubic propagator G{p) = l/\p\^ in real space. To this end, we introduce an 
infrared cutoff fi as follows 



(fp 



1 



'|p|>M (27r)3 |p|3 27r2 
where ci(A) is the cosine integral: 



sin(/i|x|) 



ci(/i|x|) 



(52) 



ci(A) 



°° COSf 



A V 



dv. 



(53) 
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As // — > we have 



[l-7-ln(/x|x|)] + 0(/i), 



(54) 



where 7 is the Euler-Mascheroni constant. For x = 0, on the other hand, Gn{x) 
is ultraviolet divergent and becomes 



GM(0) = ^ln (^) + const + 



(55) 



where A is an ultraviolet cutoff. 



Let us consider now the correlation function 



\j=i 



- J{x)ip{x)]}, 



(56) 



where J{x) = i J2j Qj^i^ — ^j) and Zq is the above functional integral for J = 0. 
Integrating out (f, we obtain 



exp 



\j=i / 
Using Eqs. (54) and (55), we obtain 



]-Y,G^{xi-Xj)qiqj 
^ i,3 



(57) 



^Gfj,{xi -Xj)qiqj^ 



27r2 



{^q)j (ln/x + 7-1) 



- In A + ^ qiqj In \xi - Xj 

i i^j 



(58) 



Thus, as /i — > the only nonzero contributions to (57) satisfy the neutrality 
condition for the charge Y^iQi = 0- The expansion in Eq. (58) is essentially the 
same as in the d = 2 case, except for the l/27r^ factor instead of a l/27r, and 
minor differences in the constants. 

The ultraviolet divergence of the phase field Ui{x) = e^nf^^'' is removed by 
introducing a wave function renormalization Q such that 



Ui{x) = Cl^'^Ui^R{x), 



(59) 
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with Ui^R being the renormahzed counterpart of Ui and 



'A' 



(60) 



Therefore, if we speciahze to the case where qi — ±|g| for all i, the renormahzed 
two-point correlation function is given by 



It follows that the dimension of Ui is just g^/(47r^). 



(61) 



Due to the above analysis it is now easy to see how Zq renormalizes in the 
ASG model. Note that the model is super- renormalizable, just as the ordinary 
sine-Gordon model. Thus, the renormalization of Zq is achieved by taking into 
account only tadpole contractions of cos (p. We obtain 



where 



(62) 



Furthermore, we have the RG function 



(63) 



The renormahzed n-point correlation function ^'^"^ satisfies the following Callan- 
Symanzik equation 



1 d \ 

Dimensional analysis, on the other hand, gives 



(65) 



g^''\p,,t,z)^o, 



(66) 



where 3(1 — n) represents the mass dimension of Using Eq. (66) in (65), 
we obtain 



d , . n f I \ 
Pi-QT, + 3(^ - 1) - + ~ 2^^J ^ 



d_ 
dz 



g^^\pi,t,z)^o. 



(67) 
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For j9 = we have 



(6 - v^)z ^ ^ = [6(1 -n) + nr/^]6;(")(0, t, (6J 

which gives the following scaling relation for small z 



(0, t, Z) ~ ^[6(l-«)+"^^l/(6-r?^) . (69) 

Also, it is clear from Eq. (67) that the scaling behavior of the mass scale is 



~ ^2/(6-r,^)_ (7Q) 

The momentum space behavior of is ~ Xjp^-'^v and therefore Q'^'^^ becomes 
singular in the ultraviolet if 3 — r^^ < 0. This happens for t = Gtc"^. For t = Gn"^ 
the mass scale behaves like z'^^^. This is an important difference between the 
usual sine-Gordon model in two dimensions and the CPSG model in three 
dimensions. The mass scale in the usual two-dimensional sine-Gordon theory 
behaves linearly in z when the singular short-distance behavior is reached. 
There, this behavior is important for the fermionization of the model, which 
establishes the equivalence between the sine-Gordon model and the Thirring 
model in two dimensions l4C 



From Eqs. (69) and (70) we see that Q^'^\0,t, z) and vanish for t = tc = 
127r^. The interpretation of this result closely parallels the one in the usual 
sine-Gordon model. For instance, it tells us that at t = t^, the operator cos<^ 
is marginal, and means that further renormalizations are necessary at t = t^. 
The situation exactly parallels the two-dimensional case where a thorough 
analysis was carried out by Amit et al. For t > tc the anomalous sine- 
Gordon model Eq. (51) is no longer renormalizable. These results follow from 
the observation that the dimension of the operator cost/? is just Thus, 
/ (i^xcos(/5 has dimension ?7^/2 — 3, which means that z has an effective di- 
mension of 3 — ?7<^/2. Therefore, the interaction is relevant for t]^ < Q or t < 
thus generating a mass. It is marginal if r^,^ = 6 and irrelevant for 77^ > 6, or 
t > tc, meaning that the theory is massless. Hence, there is a phase where the 
field has a mass and another one where it is massless, implying the existence 
of a genuine phase transition in the model Eq. (51). This follows from the fact 
that a mass changing from a finite value to zero on a finite interval of coupling 
constants must do so in a non-analytic fashion. This conclusion is one of the 
main results of this paper. Note, however, that since the above discussion is 
basically a spin wave analysis and suffices to show that a phase transition 
exists, it does not elucidate the character of the phase transition. In order to 
understand the phase transition we have to account for the topological defects 
in the theory 42|, and this is the purpose of the next subsection. 
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4-2 Kosterlitz-Thouless-like recursion relations for the anomalous sine-Gordon 
model 



The above discussion strongly suggests the existence of a phase-transition 
in the model defined by Eq. (51). However, as we have already mentioned, 
the cosine interaction becomes marginal at t = tc- This means that it is not 
true that P{t) = for all values of t. The analysis of the previous subsection is 
neglecting the monopole fiuctuations which would lead to a renormalization of 
t. This situation is well known for the lo gari thmic interaction in two dimensions 
and leads to the KT recursion relations 4^ Similar arguments can be used 



m our case. 



Let us define the dimensionless coupling y = z/fj,^. Using Eq. (62) we obtain 
the fiow equation 



The above equation can be derived in another way, which is useful for the 
purposes of this subsection. Let us consider again the monopole action in Eq. 
(49). We can write the partition function of the monopoles as 



^mon ^ ^ exp 
{n{x)} 



/ d^x / d^x'n{x)G{x - x')n{x') 

J\xi\>a J\x'^\>a 



(72) 



where a is a short distance cutoff. Using Eqs. (54) and (55), we rewrite the 
above in the following form 



{n(x)} 



Xi\>a J\x^\>a 



-2ti''K / (fx / i{'x'n{x)G(x - x')n{x') 



J \xi\>a 



(73) 



where 



'5W = -2^'°T- (") 

The prime on the summation sign in Eq. (73) indicates that the charge neu- 
trality constraint implied by the large distance limit is enforced. If we assume 
small yo such that configurations with zero or one pair of monopoles are dom- 
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inant, we obtain 



7 ~ 

•^mon ~ 



1+1/0 / 



d X 



\x[\>a \X 



1 



(75) 



i\>a J\x'^\>a \X — X ^ 

If we change the short distance cutoff in the integrals as a ^ ab, we see that 
the form of Eq. (75) is unchanged provided x and x' are rescaled in such a 
way as to restore the previous integration region and yo is changed according 



y = yob' 



3-K 



(76) 



If we define / = In 6, we obtain 



dy 
dl 



{3-K)y. 



(77) 



Recalling that t = An'^K, we see that Eq. (77) is precisely Eq. (71), except 
for the sign, which is due to differences in the cutoff procedure. Eq. (77) is 
analogous to the corresponding flow equation for the fugacity in the ordinary 
KT transition 0; 4^. In that case we find instead dy/dl = (2 — iTK)y. The 
factor 2 in the usual KT case reflects the dimensionality. In our case we have 
a factor 3 instead (and also just K rather than irK). 

It is also possible to derive recursion relations involving the fugacity of the 
problem in arbitrary dimensions to lowest orders in the fugacity for the d- 
dimensional Coulomb gas with a power law interaction 



V{x) 



(4vr)rf/2 




2-d 



1 



(7J 



This problem was considered by Kosterlitz 3| , who also obtained the flow of 
the stiffness. The result is 



dK- 



dl 



{2-d) K-\ 



(79) 



where f{d) = {d 
flow equations. 



dy 
dl 



d-2n^f{d)K 



2)T[{d - 2)/2]/(47r)^/2_ Por = 2 this reduces to the KT 



However, we see that for d = 3 the recursion relations (79) and (80) do not 
have a flxed point and therefore no phase transition happens in this case. In 
case of Eq. (49), on the other hand, we have an anomalous Coulomb gas whose 
potential is logarithmic in three dimensions. It is thus plausible to conjecture 
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that we would have a flow equation for the stiffness similar to the d = 2 KT 
case. As we shall see, this is indeed the case. 

For d = 3, Eq. (80) coincides with our Eq. (77). However, since the potentials 
are different we should in fact not expect the same recursion relation for the 
fugacity. This suggests that the "spin wave" picture of Section IV-A is not 
giving the correct flow for the fugacity. Note that for d = 2 discussed in jij], 
the spin wave analysis does in fact give the correct flow for the fugacity{see 
Appendix B) . For a logarithmic potential in ci = 3 there are some subtleties. 

Let us consider a problem with a potential like 



V{x) 




2-d+riA 



- 1 



Here we have taken into account the effect of anomalous scaling due to matter- 
field fiuctuations in our original problem. A logarithmic interaction corre- 
sponds to the case d = 3 and r]A = 1, which is the case which eventually will 
be relevant for us. Strictly speaking, the duality scenario in Section HI is valid 
only at d = 3. However, as far as the scaling behavior is concerned, it is use- 
ful to continue to the whole dimension interval (2,4), while keeping the same 
e-tensors. This dimensional continuation procedure is reminiscent of the one 
considered in some RG studies of Chern-Simons theories The recursion 
relations we obtain are given by (see Appendix B) 



dK-^ 

dl 
dy 

dl 



d-riy- 27i^f{d)K 



(82) 



where r]y is the anomalous dimension of the fugacity which is given by 



Vy 



Va 
2 



i-d 



(83) 



and 



^^"^^ ~ 2v^{47iY/^T{1 + 7]a/2) ■ ^^^^ 
Hence, for the case of a logarithmic interaction in three dimensions, which 
corresponds to ?7a = 1, the recursion relations for the fugacity and the stiffness 
have a similar structure as the standard Kosterlitz-Thouless recursion relations 
one obtains in the two-dimensional case 4^ 4^ . The main difference is in the 
recursion relation for the fugacity, which has an anomalous dimension rjy = 
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1/2. Note that the second term in the equation for K~^{1), which prevents fixed 
points of the equations (82) from being obtained, is absent for a logarithmic 
potential in any dimension. 

When rjA = 0, which corresponds to neglecting the effect of matter fields in 
the original gauge theory, we have riy = 0. Our recursion relations then reduce 
to the ones given in Eqs. (79) and (80) obtained in by a very different 
method than we employ in Appendix B. Moreover, we have also derived Eqs. 
(82) along a different route than that used in Appendix B, namely by the 
method employed in Q]. This constitutes an important consistency check on 
our calculations. For the case where r]A = 0, the absence of a phase transition 
refiects the permanent confinement of electric test charges in the usual three- 
dimensional compact QED j^. 

We see that the fiow equation for the fugacity obtained in Eq. (82) does 
not agree with the result of our "spin- wave" theory, which leads to Eq.(71) 
or, equivalently, Eq. (77). The reason for this is that an anomalous scal- 
ing dimension rjy for the fugacity is induced by the renormalization of the 
stiffness. Indeed, in Appendix B we show that a potential like (78) leads 
to an additional scaling transformation in the effective stiffness of the form 
K{1) e'^'^~'^'^^^^^K{l). If rjA 7^ 0, this is compensated in the effective fugacity 
by the scaling transformation, y{l) e~^y^y{l). In the case of the Coulomb 
gas, where rjA = 0, the spin-wave analysis gives the right answer, Eq. (80), 
as can easily be seen by working out a Callan-Symanzik RG analysis in the 
sine-Gordon theory (15) for arbitrary dimensions. Thus, deviations from an 
ordinary type of Coulomb potential in rf-dimensions lead to an anomalous 
dimension to the fugacity, Eq. (83), which cannot be obtained by spin- wave 
theory. 

The important point to note here is that a fixed point of the recursion relations 
Eqs. (82) for c? = 3 exists for the stiffness and fugacity in the limit of zero 
fugacity, so the problem scales to the weak coupling limit. Hence, the problem 
is selfconsistently found to be amenable to a KT-type of phenomenological RG 
analysis. It is not necessary to calculate to higher order in y to determine the 
fixed point. This demonstrates that the phase transition established above 
is of the KT type. This has some resemblance with the results of a rather 
remarkable paper by Amit et al. [4^, which also finds a KT transition in 
a three-dimensional Coulomb gas with logarithmic interaction between point 
charges (see Appendix A). In their case, the logarithmic interaction between 
the point charges in three dimensions did not have its origin in anomalous 
scaling dynamically generated by matter-field fiuctuations, but originated in 
anisotropic higher order derivative terms in an underlying field theory that 
were put in by hand. This anisotropy ultimately induces a dimensional reduc- 
tion. 
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In four dimensions, we have rjA = and extrapolating the above results it is 
clear that no fixed points of the above recursion relations can be found. Indeed, 
the above analysis no longer applies and no KT topological phase-transition 
occurs. This is so because by dualizing a compact Maxwell Lagrangian in four 
dimensions, we obtain a non-compact Abelian Higgs model jSSf , which cannot 
be brought onto the form of a Coulomb gas. The transition in this case is 
known to be of more conventional second or first order type j^. 

Finally, we note that in three dimensions there is a universal jump in the 
stiffness parameter at the transition, analogous to what is known in the 2d 
case (4^. In units of Eqs. (82), this jump is determined by dimensionality and 
the anomalous scaling of the fugacity. 



5 RG functions of Abelian Higgs model and KT transition 

In this Section we show how the RG functions and fixed points in the Abelian 
Higgs model are related to the KT-like transition described in the previous 
Section. In particular, we shall use the critical coupling tc to fix an a priori 
arbitrary constant that enters into the computation of the critical exponents 
for the Abelian Higgs model. This in our view improves on a scheme previously 
used (2^, where a corresponding constant was fixed by appealing to numerical 
results for the value of the Ginzburg-Landau parameter k which separates 
first- from second-order behavior ^ . In our approach, the parameter (denoted 
r below) is fixed from our theory of the critical behavior of the compact case, 
which we have argued in the Introduction to be the same as for the non- 
compact Abelian Higgs model at infinite bare gauge coupling. Before doing 
this, however, a few preliminary remarks are in order. 

The Abelian Higgs model is manifestly a two-scale theory. Indeed, the gauge 
field becomes massive due to the Higgs mechanism. Thus, in the ordered phase 
we are left with two mass scales, the Higgs mass m and the gauge field mass 
m^. From these two mass scales we obtain the Ginzburg parameter k = m/niA. 

^ In an early Monte Carlo simulation, a tricritical value Ktri = 0.4/\/2 was found, 
14^ 1 . This is the value used in the ad hoc scheme of Ref . . More recently, a large- 
scale Monte Carlo simulation improved on this value, finding Attn = (0.76ib0.04) /\/2, 
ji^ . This is in sur pris ingly good agreement with an early analytical result Ktri = 
0.798/\/2, see Ref. |27l |. Using this improved value for Ktri iu the /^-functions of Ref. 
[2]| . the critical exponent u obtained would be = 0.53. This is quite far from 
the correct 2>DXY value vxy = 0.67, as well as from the ?>DXY one-loop value 



Kr = lim K{1) 



d-r]y 
27r2/(rf) 




u = 0.625. 
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Due to the existence of two mass scales in the problem, we have very distinct 
situations depending on whether k -C 1 or k ^ 1. For k <^ 1 vortex lines, 
which are the topological defects of the matter field, attract each other. This 
corresponds to a type I regime, while for k ^ 1 we have repulsive forces 
between vortex lines, which corresponds to the type II regime. This two-scale 
behavior survives in the disordered phase, though in this case ttla = 0. 



We shall consider the calculation of RG functions for the massless theory, but 
using two renormalization scales 
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In order to see the influence of the two 
mass scales appearing in the ordered phase, on the massless theory, we define 
the dimensionful couplings at different renormalization points, u at fi and 
at fi. Let us define the ratio r = n/jl. By rewriting e^(/i) in terms of /z, we 
obtain the one-loop /3-functions for any fixed dimension d G (2, 4] and an order 
parameter with N/2 complex components [50] 



/3„ = (4 - d)[-a + rNA{d)a\ 



(86) 



(5g = {A-d)^^-g + B{d) 



-2{d- l)ag + 



57) 



where 



A{d) 



r(l - d/2)T^{d/2) 
(47r)'^/2r(rf) ' 



(88) 



B{d) 



r(2-rf/2)r2(d/2- 1) 



9) 



(47r)^/2r(ci-2) ■ 

From Eq. (86) we see that 7^ = r{4: — d)N A{d)a. By considering d = A — e and 
expanding for small e, we recover the well known e-expansion result [51] if we 
take r = 1. In our fixed dimension approach r is an arbitrary parameter that 

When d = 3 and iV = 2 



21 



is usually fixed by imposing additional conditions 
we have the fixed point a*(r) = 16/r. In the context of the compact Abelian 
Higgs model we fix the value of r by demanding that K^. should correspond 
to a r = Tc, with K = l/a* at one-loop. If we use the spin-wave estimate 
Kc = 3 (which corresponds to = 127r^), we obtain then that = 48 and 
thus a* = 1/3. On the other hand, if we use the estimate from our KT-like 
recursion relations, we have Kc = 5/2 and therefore a* = 2/5. In order to 
check the quality of these matchings, we compute the critical exponents of the 
three-dimensional Abelian Higgs model in d = 3. The critical exponent z/ is 
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given by the fixed point value of tlie RG function 



where 

d\xiZ, 

with Zm being the mass renormahzation and 



7m = f^—Qj^ 70, (91) 



7* = .-^. (92) 

At the fixed point 7^ gives the value of the critical exponent 77. At one- loop 
order, we have 



7m = ^^j^, 70 = ~- (93) 

When K(. = 3, the fixed point for the coupling g which corresponds to infrared 
stability is given by = 2(7 + 2v/TT)/15. Therefore, we obtain u ^ 0.615 and 
r] = -1/12. Using = 5/2, we obtain g, = 4(6 + y3T)/25. The critical 
exponents in this case are u ~ 0.61 and 1] = —1/10. Both estimates are close 
to the one-loop value of the XY model, uxy ~ 0.625. From duality arguments 
we expect indeed a XY value for the exponent u [s^ . 



6 Summary and Discussion 



In this paper, we have considered the Abelian Higgs model in 2 + 1 dimensions 
both for the non-compact and compact cases, with matter fields in the funda- 
mental representation. We have performed a duality lattice transformation on 
these models, emphasizing the features that set them apart as well as those 
they have in common. A major difference lies in the fact that in the dual for- 
mulation, the non-compact case has stringent constraints V ■ Ij = imposed 
on the topological currents of the system, while in the compact case V ■ li can 
take any integer value, i.e. the currents are unconstrained for the case where 
the matter field is in the fundamental representation. This effectively makes 
the dual non-compact case a much more strongly interacting system of topo- 
logical currents, and this is why phase transitions are more easily brought out 
compared to the compact case. As a result, we have seen that there is one limit 
of the LAH model where the non-compact case exhibits the IXY transition, 
while the compact case is an exactly soluble discrete Gaussian model with 
apparently no phase transition. 



28 



A major part of the paper (Section III and IV) has been devoted to estab- 
hshing that, despite the absence of any phase transitions with a local order 
parameter in the compact topological phase transition nevertheless is 

found in the interior of the phase diagram of the model. A key ingredient 
is the renormalization of the gauge-field propagator of the problem due to 
critical matter field fluctuations, Eq. (46). With no matter fields present, the 
topological defects of the gauge field, which are monopole configurations, in- 
teract with a 1/R potential in d = 3. In the presence of matter fields, taking 
into account their critical fluctuations, the resultant effective gauge theory 
may be described as an overall neutral plasma of charges that interact with 
a logarithmic potential in c? = 3, Eq. (49). A field-theoretical formulation of 
the action given in Eq. (49) yields an anomalous sine-Gordon (ASG) model, 
Eq. (51). A renormalization group analysis of this model based on the Callan- 
Symanzik equations shows that the theory is massive below a critical value of 
the couphng constant. This by itself suffices to conclude that a phase transi- 
tion exists. We then go on to show that the problem is amenable to an analysis 
based on KT-like recursion relations, Eqs. (82), derived for a d-dimensional 
gas of point charges interacting with a pair-potential which in a certain limit is 
logarithmic. In this limit, the recursion relations we derive for the stiffness and 
fugacity of the problem reduce to equations which are similar in structure to 
the well known Kosterlitz-Thouless recursion relations obtained for the two- 
dimensional Coulomb gas, but with a modified equation for the fugacity due 
to an induced anomalous scaling of it. This anomalous scaling in the fugacity 
accounts for deviations from the ordinary Coulomb gas case in d dimensions. 
The change in the equation for the fugacity shows that the stiffness and the 
fugacity of the problem mutually influence each other under renormalization 
in a manner which is different from the case of a logarithmic pair-interaction 
in (i = 2. As a consequence of this, the universal jump in the stiffness at the 
transition is then given, in appropriate units, by the dimensionality of the 
system and the anomalous scahng of the fugacity, Eq. (85). 

In Section V, we have seen that the deconfinement phase transition we find in 
the compact case, with a critical coupling tc, allows us to fix a parameter ap- 
pearing in the evaluation of the critical exponents of the non-compact Abelian 
Higgs model. This represents an improvement on previous schemes to fix this 
parameter. 

We close with a few remarks on unsolved problems. When only fermionic fields 
are coupled to the massless gauge field (spinor QED^), then we again obtain 
a P function for the renormalized gauge coupling as given in Eq. (4), but 7^ 
in the equation now depends only on one coupling constant, a, not two as 
in the bosonic case. Then we do not have the freedom to tune parameters 
of the model to drive it through a phase transition of the type described in 
Section IV. The analysis of Section IV may be carried through as before, but 
the point is that the fixed point coupling, a — does not depend on any 
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second coupling constant g, this simply does not appear in the theory. Instead, 
a* depends on the number of fermion flavours only. In principle there thus 
exists a critical value N = where the compact version of the model with 
fermionic matter, also goes through a deconfinement transition. The confining 
phase corresponds to < A'^c- It is highly controversial what this critical 
value is. A simple one-loop renormalization group calculation gives Nc = 24 
(20I 1 in agreement with an earlier result by loffe and Larkin obtained by a quite 
different method However, we may in fact expect that the actual value is 
much smaller than this. Marston has calculated the same number using one- 
instanton action and finds = 0.9 [53j. The important point here is that 
whatever the precise value of Nc is, the interaction between the monopoles is 
always logarithmic. 



Also, in the fermionic case there is a subtlety in that another type of instabil- 
ity, absent in the bosonic case, could intervene to destroy the deconfinement 
transition. Fermions can in principle undergo a spontaneous chiral symmetry 
breaking (S^SB) js^]- This happens when the number of fermion flavours is 
less than some critical value, A"ch say. This means that a fermion mass is dy- 
namically generated for A^ < A"ch- The precise value of A'ch is presently also 
a matter of debate. One estimate from the Schwinger-Dyson equation gives 
A'ch = 32/7r^ 55|. This result is confirmed by Monte Carlo simulations finding 
A'ch ~ 3.5 [56|. Another analytic calculation gives A'ch = 128/37r^ ji^l- A re- 
cent estimate based on a new constraint on strongly interacting systems gives 
-^ch < 3/2 5^1 . This is quite consistent with the most recent numerical results 
we are aware of 5^|, where no signs of S^SB is found for N > 2. Thus, there is 
no consensus on the precise value of A^ch- The calculation of Nc assumes that 
the fermions are massless. Thus, if Nc = 24 as in Refs. ji^ and jl^l, then a 
deconfinement transition will take place since the fermion mass is generated 
at a much lower value of A^. With massive fermions present our anomalous 
three-dimensional compact QED scenario does not apply because the Maxwell 
term does not become irrelevant anymore. In such a situation the results of 
Polyakov apply and there is permanent confinement of electric test charges. 
This would be the case for the value Nc = 0.9 obtained by Marston [s^, which 
lies below all estimates of A'ch- In this case the deconfinement transition does 
not happen. 



Physically, S^SB in spinor QED^ has important consequences in the physics 
of high-Tc cuprates. As we mentioned in the Introduction, spinor QED^ with 
a compact gauge field emerges as a possible low energy description of the fluc- 
tuations around the flux phase in the quantum Heisenberg antiferromagnet 
0. In this context, the djTiamical mass generation is associated with the spin 
density wave (SDW) instability. Thus, gauge field fluctuations could in prin- 
ciple restore the Neel state. The physical number of fermion components in 
this case is A^ = 2. Spinor QED^ also emerges by considering the low energy 
physics of the d-wave superconducting state in the pseudogap phase of the 
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high-Tc cuprates 6^1 . In this case, however, the gauge field is non-compact 
and there is an inherent anisotropy in the Lagrangian. There also, S^SB is 
responsible for the onset of SDW as half- filling is approached ^]. The phys- 
ical number of fermion components in this case is again N = 2. Therefore, 
in these theories it is essential that A^^ch > 2. If the most recent estimate for 
A^ch is correct [s^, this could have serious implications for the validity of the 
different spinor QED3 scenarios discussed above. In the case of the spinor 
QED^ description of the pseudogap phase, the inherent anisotropy could pos- 
sibly affect the value of N^t- However, results presented thus far indicate that 
at least weak anisotropy will not affect A^^ch obtained in the isotropic case [sj]. 
Moreover, when studying effective theories of undoped high-Tc cuprates, we 
have argued in the Introduction that the relevant theory to study is fermions 
coupled to compact U{1) gauge- fields. Hence ^ it is of importance to revisit the 
problem of how monopoles affects S^SB |63|. Finally, we note that a recent 
provocative paper by Wen 0] states that there exists a principle of quantum 
order which may prevent fermions from dynamically acquiring a mass even in 
the presence of strong coupling to gauge fields. Hence, it seems to us that a 
renewed effort in numerical computations of Nc\i in 2 + 1-dimensional gauge 
theories coupled to fermionic matter, including the effects of compactness and 
anisotropy, would be very timely. 
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A KT-like transition in three dimensions in an anisotropic sine- 
Gordon theory 



While considering a class of globally symmetric self-dual Zp^ models in the 
A ^ 00 limit, Amit et al. |46[ arrived at the following anisotropic three- 
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dimensional sine-Gordon action containing higher derivatives: 



ANISG 



(A.l) 



_2ti II 2t2 

where 9y = (9^ + 9^. As pointed out in Ref. ji^, the above model has a KT 
transition in three dimensions. Indeed, it is easy to see that the propagator is 
logarithmic at large distances. Note, however, that anisotropy and the higher 
order derivatives in the parallel direction are essential, and the system effec- 
tively shows two-dimensional behavior by dimensional reduction. This is in 
contrast with our genuinely three-dimensional KT-like scenario. 



B KT-like recursion relations 



In this appendix we derive to lowest order in the fugacity the recursion rela- 
tions for the scale-dependent stiffness parameter K{1) and fugacity y{l) given 
in Eqs. 82 for a d-dimensional plasma where the bare pair-potential is given by 
Eq. (81), which reduces to a logarithmic potential when d = 3. The starting 
point will be a low-density approximation for a dielectric constant of this sys- 
tem. We closely follow a method for doing this introduced in . Introducing 
the solid angle in d dimensions Vld = 2n^^'^ /T(d/ 2) and the density of dipoles 
in the fluid by n^^, a low-density approximation for the dielectric constant is 
given by 

e = l+ndnda (B.l) 

where a here denotes the polarizability of the medium, a standard linear- 
response analysis gives a = Att'^K < > /d and < > is the mean square of 
the dipole moment in the system. To compute this, we need the low-density 
limit of the pair-distribution function n^(r) of the plasma, which is readily 
obtained from the grand canonical partition function S expanded to second 
order in the bare fugacity and replacing the thermal de Broglie wavelength 
by a short-distance cutoff ro, as follows 



nHr) = ^ e-^-^^^. (B.2) 



In this way, we may now go on to express a scale-dependent dielectric constant 
as follows 

e(r) = 1 + r ds n^is). (B.3) 

Note however, that in Eq. (B.3), a mean-field approximation is understood to 
be used by replacing the bare potential V in n^(r) by an effective potential 
U{r). This effective screened potential must be selfconsistently determined by 
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demanding that it gives rise to an electric field in the problem given by 



?1 . E^r) ^ (B.4) 

or £[r) r p 

where p = 2 — d + r]A and f{d) is defined in Eq. (84). Such a mean-field 
procedure has been consistently used with success in the 2d case, and the origin 
of the success lies in the long range of the In-interaction. In higher dimensions, 
such a procedure will work even better since the logarithmic potential is felt 
over even longer distances due to extra volume factors. 

Let us introduce a logarithmic length scale I = ln(r/ro) along with the new 
variables 



5(roexp/) 

x{l)=A7i^KU{roexpl) (B.5) 

Here, x{l) is determined selfconsistently by integrating the effective field E{r). 
Then we get from Eqs. (B.3) and (B.4) 

r{l) = r(0) + t dv e('^+2).-x(.)^ (g g^ 

dvQ Jo 

and 

x{l) = x{0) + f{d) f dv (B.7) 
Jo t{v) 

From Eqs. (B.6) and (B.7), we may derive coupled renormalization group 
equations for r(/) and x{l). However, in order to obtain equations that have a 
form more similar to equations that have appeared in the literature on the d- 
dimensional Coulomb gas ji^l, we introduce a new variable K{1) representing 
a scale dependent stiffness constant, as follows 

Thus, we see that the effect of a nonzero p on the stiffness amounts to a scaling 
change K{1) — eP^K{l). Using Eq. (B.7), we have that 

^ = 4vr7(c^)i^(/). (B.9) 
Differentiating K~^{1) with respect to / and using Eq. (B.6), we obtain 

= -P K-\l) + e[('^+-^)-^(0] (B.lO) 
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Prom this expression, we define a scale dependent fugacity y{l) given by 



y^'d C e^(d+2-p)l-x{l)]/2 



2/(0 = r-,Sd-2+p)/2 ■ 



Thus, we see exphcitly that the renormahzation of K{1) in principle influences 
the flow equation for y{l), which is obtained by differentiating with respect to 
/ and using Eq. (B.9) 



dyjl) 
dl 



d-r]y-2n'f{d)K{l) y{l), (B.12) 



where rjy = {d — 2 + p)/2. Eqs. (B.IO) and (B.12) are precisely Eqs. (82). 
On the other hand, the Callan-Symanzik approach of Section IV-A, which 
basically ignores the influence of the renormalization of K{1) on the structure 
of the flow equation for y{l), yields as we have seen Eq. (77). We have already 
remarked in Section IV-B that this type of approach gives the correct answer 
only if there are no deviations from the Coulomb potential case, that is, we 
need p — 2 — d. Note that in the usual KT transition we would have p — rjy — 0. 



C Screened effective potential 



In this appendix, we derive the asymptotic long-distance behavior of the 
screened effective interaction U{r) introduced in Appendix B, for the case 
p = 0, corresponding to o? = 3 and rjA = ^- We start from the recursion 
relations, written on the form 



dK 



-1 



dy_ 
dl 



-m 



y- 



(C.l) 



From Eq. (B.8), wc have that K ^(l) 
the variable T(/) defined by 



r(/) in this case. Next, we introduce 



5r(0/2-l 

5r(/)/2 

«^r(0-l, (C.2) 

where the latter approximation is asymptotically exact close enough to the 
transition. In terms of this, the flow equation for the fugacity may be written 
on the form 



dy'jl) 
dl 



5 T{1) y\l), (C.3) 
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On the other hand, we have 



dT\l)^,^^ dr(l) 
dl dl 

= 5 T{1) y\l), (C.4) 



and hence we have 



y\l)-T'(l)^±u;', (C.5) 

where lv is some positive number. We are interested in the quantity hm^^oo 
for the case where y'^{l) — T'^{1) < 0, and T{1) < 0, this wiU be the regime 
where the fugacity scales to zero. In this case we choose the negative sign on 
the r.h.s. in Eq. (C.5). From the flow equation for K~^(l) we flnd 



^^(0-5 . 5 
dl -2^^^'' 2 



(C.6) 



This is solved to obtain, introducing u = {5/2) lu I + 9, 



T{l)^-uj cothw 

= (C.7) 
smhw 

where uj and 9 are integration constants that are uniquely determined from 
the initial conditions on r{l) and yil), i.e. by the bare coupling constants of 
the problem as follows 



y^(0)-r2(0) = -a;^ 

r(o) 



cosh^. (C.^ 



y(o) 

Prom the expression for T(/), using Eq. (C.2), we obtain 



2 

t{1) ^-{i-oj cothw). (C.9) 
5 



Since t{1) > 0, this puts restrictions on the constants u) and 9, and the most 
severe limitations on cu in terms of 9 is given by 

l-a;coth^>0. (C.IO) 

Using Eq. (B.9) and the fact that K{1) = 1/t(Z), we have 

dx{l) 5/2 



dl 1 — a; coth u 



(C.ll) 
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From Eq. (CIO), we see that dx{l)/dl > 0. This is an important result, since 
it immediately reveals that, in the regime y'^{l) — T'^{1) < we consider here, 
the logarithmic bare potential V{r) cannot possibly be screened into a power 
law potential with a > 0, since in that case we would have dx{l)/dl < 0. 
However, for all / we have 



d^x{l) 



5 uo/2 



< 0. 



(C.12) 



dP ysinh-u — C(j cosh-u^ 

Introducing cu^ = 1 ± cj, Eq. (C.ll) is straightforwardly integrated to yield 



x{l) -a;(0) 



1 



-cj+Z + In 



-2u 



+ UJ- 



(C.13) 



From this, it follows that for r ^ vq the effective potential behaves asymptot- 
ically as 

f/(r) ~ln(r/ro). (C.14) 



D Exact equation of state for the rf-dimensional In-plasma 



The equation of state for a d-dimensional In-plasma with no short- distance 
cutoff, may be obtained via a simple scaling argument, previously applied to 
the two-dimensional case The configurational integral in the canonical 



partition function is given by 

Q = /^■■■J^ d'^vi ■ ■ ■ d^r^N exp[t ^ ln(r,,)] (D.r 



where qi = ±1, and we assumed that we have 2N particles in the system, N 
with charge qi = 1 and with charge = —1, ^2^=1 Qi = 0. Here, V = L'^ is 
the volume of the system. Introduce new dimensionless variables Rij = Vij/L 
where r^j = | r j — r j | , in which case the configurational integral is given by 



Q = L^^'i f ... f d'^Ri ■ ■ ■ rf^i?27V exp(t V qiqj HRijL)) 
Jo Jo 

= L'^''exp[iJ2q,q,ln{L)] I, (D.2) 

i<j 

where the integral / is independent of volume. Now note that 



2 

i<j 



■■-2N 



2N 



(D.3) 
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Then we obtain 



Q _ jJ2Nd^-iNln{L) j _ jJ2Nd-iN j 

^y2N-iN/d j_ P_4) 

From this, we obtain the equation of state involving the pressure 

f N 

tpV = 2N -. (D.5) 

a 

Note that the pressure vanishes when i = = 2d. A prerequisite for the 
vahdity of the above analysis is that the quantity / must be finite, otherwise 
the scahng of variables that lead to the equation of state is meaningless. In 
fact, / is not always finite. Consider again the integrand in Q, which is given 
by a product of factors 

Any factor with = —qj will be singular when Vij = 0, which is possible 
in the absence of a short-distance cutoff. To investigate whether or not this 
singularity is integrable, consider the integral 

J drr"^-^ r"* (D.7) 

This is finite only if 

d-i>0 (D.8) 

This means that the equation of state Eq. (D.5) makes sense for t < tc = d, 
note that for all dimensions d, to = Ste- 
in two dimensions, it is known that the negativity of the pressure occurs 
at a temperature that coincides with the KT vortex-antivortex unbinding 
temperature, and that there is a phase transition at twice this temperature. 
It is amusing to note here that in the three-dimensional case, the pressure 
vanishes at tc — 127r^, after having reintroduced t = t/An'^. This is precisely the 
critical coupling we found in Section IVA from the Callan-Symanzik equations. 
In addition there is again a phase transition at precisely half the value of 
this coupling constant, where the pressure becomes that of an ideal gas of N 
particles. In arbitrary dimensions, this persists, the phase transition to an ideal 
gas of A'" particles always happens at half of the value at which the pressure 
vanishes. This phase transition, which is a collapse of an overall charge- neutral 
plasma of A^ = +1 charges and A^ = — 1 charges into an ideal gas of A^ 
particles, occurs because of the lack of a short- distance cutoff in the system 
we consider in this appendix. 
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E Duality in the Abelian compact Higgs model with a Chern- 
Simons term 



For completeness, we present in this Appendix the duality transformation of 
the LAH with a Chern-Simons term added 3|. Compact gauge theories with 



11 



Chern-Simons term added are relevant in studies of chiral spin liquid states 
when spinor states have been integrated out. Such theories have been argued 
to exhibit a deconfinement transition 67; [g^. The compact LAH mode, i.e. 
Aifj^ G (— vr, vr), with a Chern-Simons term has the action 



cs ■ 



E 



•7, 



(E.l) 



Let us introduce auxiliary fields a^, bj, Aju, and Uiu, such that 



P 2 

2^'^ 2e2"^ ■ '2 
+ icTifii^fiux'^uAix - 2nNi^ - bi^)] . 



(E.2) 



Next we introduce integer valued fields rrij^ and Mj^ via the Poisson formula: 



i 



(E.3) 



Integration of 9i and Ai^ give the constraints enforced by delta of Kronecker 



V ■ m,; = 0, 



(E.4) 



Summing over m^ gives 



V X M, = m, . 



(E.5) 



Qlll 



E 



^a^ + -^h] + i^Sii ■ hi - iiy X Mi) ■ a^ - iM, ■ b^ 



(E.6) 



38 



By integrating out aj and bj we arrive at the action 



(E.7) 



where K = 4/ (7^6^ + 4/3). Using the Poisson formula to introduce a real lattice 
field hifj, and doing an appropriate rescaling of the variables we obtain finally 
the partition function 



exp[-^^f (h,,l 



(E.8) 



where 



which should be compared with Eqs. (26) and (34). Note the appearance of 
the cross-term ije^hi ■ (V x hi). When the hj are integrated out we are thus 
left with a partition of the same form as Eq. (35), but with an asymmetric 
propagator. 

If we were to consider the non-compact LAH with a Chern-Simons term added, 
and in the absence of the Maxwell term, — > oo, then this is an effective 
description of the fractional quantum Hall effect j37|;|69|. In this case we obtain 



4^(V X hi? - —hi ■ (V X hi 



+ i2n\i-h,. (E.IO) 



This is essentially the same as Eqs. (E.8) and (E.9) for the compact case (with 
no mass term for the hj-fields), but we should add an additional constraint in 
the V ■ Ij = in the partition function. 

One point worth emphasizing here, sometimes overlooked, is that the gauge- 
field hj is never a compact gauge-field, whether one starts from an original 
compact or non-compact gauge theory. In the non-com pac t Chern-Simons the- 
ory, there exists a self-dual point at a value 7 = l/27r [g^ [t^. The possibility 
of self-duality is a consequence of non-compactness, it can never arise starting 
from a compact LAH model with Chern-Simons term added. It is an intrigu- 
ing question whether the self-duality at the above particular value of 7 in the 
non-compact case corresponds to a critical point. A candidate physical inter- 
pretation of such a putative phase transition would correspond to statistical 
transmutation of the Laughlin quasiparticles of the fractional quantum Hall 
effect as magnetic field is varied, since in the context of the FQHE, the pa- 
rameter 7 depends on filling fraction, i.e. magnetic field. It is known that for 
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the half-filled lowest Landau level, the quasiparticles are fermions [71|, while 
for other filling fractions they are anyons. 
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